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I. INTRODUCTION 

B meson decays to final states containing an axial- 
vector meson (A) and a pseudoscalar meson (P) have 
been studied both theoretically and experimentally. The- 
oretical predictions for the branching fractions (BFs) of 
these decays have been calculated assuming a naive fac- 
torization hypothesis P, [1| and QCD factorization 
These decay modes are expected to occur with BFs of 
order 10~^. Branching fractions of B meson decays with 
an ai(1260) or &i(1235) meson plus a pion or a kaon in 
the final state have recently been measured [1, Q . 

The BaBAR Collaboration has measured CP-violating 
asymmetries in P° ai(1260)^7rT decays and deter- 
mined an effective value acff [6] for the phase angle a of 
the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing 
matrix [3]. In the absence of penguin (loop) contribu- 
tions in these decay modes, acff coincides with a. 

The AS = 1 decays we examine here are particularly 
sensitive to the presence of penguin amplitudes because 
their CKM couplings are larger than the corresponding 
AS* = penguin amplitudes. Thus measurements of the 
decay rates of the AS* — 1 transitions involving the same 
SU(3) flavor multiplet as ai(1260) provide constraints on 
Aa = acff — a [8|]. Similar SU(3)-based approaches have 
been proposed for the extraction of a in the tt^tt" Q, 
p=^7r=F [8], and p+p' channels [lS[iH. 

The rates of i? — > Kiatt decays, where the Kia meson 
is the SU(3) partner of ai(1260) and a nearly equal ad- 
mixture of the ii'i(1270) and _ft'i(1400) resonances p^ . 
can be derived from the rates of i? i4ri(1270)7r and 
B i^i(1400)7r decays. For B° ii:i(1400)+7r- [3 
and i?+ iiri(1400)'^7r+ decays there exist experimen- 
tal upper limits at the 90% confidence level (C.L.) of 
1.1 X 10-3 and 2.6 x 10"^, respectively [H. In the fol- 
lowing, Ki will be used to indicate both i4ri(1270) and 
i4'i(1400) mesons. 

The production of Ki mesons in B decays has been 
previously observed in the B J/tpKi, B Kij, and 
B — > Kicj) decay channels [l^. Here we present mea- 
surements of the P° — > K^Tr~ and B^ — > K^tt^ branch- 
ing fractions and estimate the weak phase a from the 
measurement of the time dependent CP asymmetries in 
B'^ ai(1260)^7r^ decays and the branching fractions 
of SU(3) related modes. 

This paper is organized as follows. In Sec. [IT] we de- 
scribe the dataset and the detector. In Sec.|Tn]we intro- 
duce the _fC— matrix formalism used for the parameteriza- 
tion of the Ki resonances. Section HVl is devoted to a dis- 
cussion of the reconstruction and selection of the B can- 
didates. In Sec.|V]we describe the maximum likelihood fit 
for the signal branching fractions and the likelihood scan 
over the parameters that characterize the production of 
the Ki system. In Sec. |Vl] we discuss the systematic 
uncertainties. In Sec. IVIII we present the experimental 
results. Finally, in Sec. Ivml we use the experimental 
results to extract bounds on |Aa|. 



II. THE BABAR DETECTOR AND DATASET 

The results presented in this paper are based on 
data collected with the EaBAR detector at the PEP- 
II asymmetric-energy e^e~ storage ring, operating at 
the SLAC National Accelerator Laboratory. At PEP-II, 
9.0 GeV electrons collide with 3.1 GeV positrons to yield 
a center-of-mass (CM) energy of ^/s — 10.58 GeV, which 
corresponds to the mass of the T{AS) resonance. The 
asymmetric energies result in a boost from the labora- 
tory to the CM frame of « 0.56. We analyze the final 
BaBAR dataset collected at the T{AS) resonance, corre- 
sponding to an integrated luminosity of 413 fb"^ and 
Ngg = (454.3 ± 5.0) X 10^ produced BB pairs. 

A detailed des crip tion of the BaBAR detector can be 
found elsewhere [l6|. Surrounding the interaction point 
is a five-layer double-sided silicon vertex tracker (SVT) 
that provides precision measurements near the collision 
point of charged particle tracks in the planes transverse 
to and along the beam direction. A 40-layer drift cham- 
ber surrounds the SVT. Both of these tracking devices 
operate in the 1.5 T magnetic field of a superconduct- 
ing solenoid to provide measurements of the momenta 
of charged particles. Charged hadron identification is 
achieved through measurements of particle energy loss in 
the tracking system and the Cherenkov angle obtained 
from a detector of internally reflected Cherenkov light. 
A CsI(Tl) electromagnetic calorimeter provides photon 
detection and electron identification. Finally, the instru- 
mented flux return (IFR) of the magnet allows discrimi- 
nation of muons from pions and detection of mesons. 
For the first 214fb-^ of data, the IFR was composed of 
a resistive plate chamber system. For the most recent 
199 fb^^ of data, a portion of the resistive plate chamber 
system has been replaced by limited streamer tubes . 

We use a GEANT4-based Monte Carlo (MC) simula- 
tion to model the response of the detector [l8|, taking 
into account the varying accelerator and detector con- 
ditions. We generate large samples of signal and back- 
ground for the modes considered in the analysis. 

III. SIGNAL MODEL 

In this analysis the signal is characterized by two 
nearby resonances, i^i(1270) and if i (1400), which have 
the same quantum numbers, I{J^) — 1/2(1+), and decay 
predominantly to the same Ktttt final state. The world's 
largest sample of Ki{\21Q) and if i (1400) events was col- 
lected by the ACCMOR Collaboration with the WA3 ex- 
periment [l^. The WA3 fixed target experiment accu- 
mulated data from the reaction K~p — > K~TT^Tr~p with 
an incident kaon energy of 63 GeV. These data were ana- 
lyzed using a two-resonance, six-channel if-matrix model 
[2(j | to describe the resonant Kmr system. We base our 
parameterization of the if i resonances produced in B de- 
cays on a model derived from the if -matrix description 
of the scattering amplitudes in Ref. [l9(. In Sec. IIII Al 
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we briefly outline the _K'-matrix formalism, which is then 
applied in Sec HITBl to fit the ACCMOR data in order to 
determine the parameters describing the diffractive pro- 
duction of Ki mesons and their decay. In Sec. IIII CI we 
explain how we use the extracted values of the decay pa- 
rameters and describe our model for Ki production in 
B {K-Kn)-K decays. 



A. 7('-matrix formalism 

Following the analysis of the ACCMOR Collabora- 
tion, the ifTTTT system is described by a i^T-matrix model 
comprising six channels, 1 = (if *(892)7r)s, 2 = pK, 
3 = ii'o*(1430)7r, 4 = /o(1370)i4:, 5 = {K* {S<d2)TT) d , 
6 = ojK . We identify each channel by the intermedi- 
ate resonance and bachelor particle, where the bachelor 
particle is the tt or if produced directly from the Ki de- 
cay. For the if*(892)7r channels the subscript refers to 
the angular momentum. 

We parameterize the production amplitude for each 
channel in the reaction K~p {K^Tr^Tr^)p as 



(1) 



where the index i (and similarly j) represents the i*^ 
channel. The elements of the diagonal phase space matrix 
p{M) for the decay chain 



ifi -> 3 -I- 4, 3-^5 + 6, 
are approximated with the form 



(2) 



2(5,,- 



2m*m4 

•n* + TO4 



(M - TO* - TO4 + iA), (3) 



where M is the Kttt: invariant mass, TO4 is the mass of 
the bachelor particle 4, and to* (A) is the pole mass (half 
width) of the intermediate resonance state 3 [21|. In Eq. 
((H), the Si parameters are offset phases with respect to 
the {K*{892)tt)s channel (^i = 0). The 6x6 if-matrix 
has the following form: 



faifaj fbifbj 

Ma - M Mb - M ' 



(4) 



where the labels a and b refer to ii'i(1400) and ifi(1270), 
respectively. The decay constants fai, fti and the K- 
matrix poles Ma and Mb are real. The pro duction vector 
P consists of a background term D [23] and a direct 
production term R 

P = (1 + tK)D + R, (5) 

where r is a constant. 

The background amplitudes are parameterized by 



D, = D 



iO 



M2 



(6) 



for all channels but (if*(892)7r)£) and usK. For the 
(ii:*(892)7r)£, channel we set = as in the ACC- 
MOR analysis [l3|. The parameters for the uiK channel 
are not fitted, as described later in this Section, and we 
set = 0. The results are not sensitive to this choice 
for the value of De- 
ll is given by 



Ri 



fpafa 



fpbfb: 



Ma-M Mb-M' 



(7) 



where fpa and fpb represent the amplitude for producing 
the states i5ri(1400) and ifi(1270), respectively, and are 
complex numbers. We assume fpa to be real. 

P-wave {£ = 1) and D-wave {£ = 2) centrifugal barrier 
factors are included in the Ki decay couplings fai and 
fbi and background amplitudes DiQ, and are given by: 



B,{M) 



q^iM)^R^ 
l + qi{MYR^ 



(8) 



where qi is the breakup momentum in channel i. Typical 
values for the interaction radius squared R^ are in the 
range 5 < i?^ < 100 GcV^^ ^ and the value R^ = 
25 GeV"^ is used. 

The physical resonances if 1 (1270) and if 1 (1400) are 
mixtures of the two SU(3) octet states Kia and if is: 

|is:i(1400)) = |is:iA)cos6l+ |iCiB)sin6i, (9) 
|is:i(1270)) = -|iCiA)sin6l+ |is:iB)cos6'. (10) 

Assuming that SU(3) violation manifests itself only in 
the mixing, we impose the following relations jl9i |: 



fai — 


^7+ cos + 1 


/ 20 1- 


sin0. 


(11) 


fbl = 


-57+ sin6' + 


Y 20 


_ cos6'. 


(12) 


fa2 = 


^7+ cos 6 — ^ 


/ 20 1- 


sin 9, 


(13) 


fb2 = 


— ^7+ sin 6* — 


V 20 1- 


- cos 9, 


(14) 



where 7+ and 7_ are the couplings of the SU(3) 
octet states to the (if * (892)77)5 and pif channels: 
((if*(892)7r)s|ifi^) = i7+ = (pif|ifiA) and 

(pif lifis) = -^^7- - - ((if*(892)7r)5|ifiB). The 

couplings for the toK channel are fixed to 1/a/3 of the 
pif couplings, as follows from the quark model [l9j . 

B. Fit to WAS data 

Only some of the if -matrix parameters extracted in 
the ACCMOR analysis have been reported in the litera- 
ture 19]. In particular, the results for most of the decay 
couplings fai and fbi are not available. The ACCMOR 
Collaboration performed a partial-wave analysis of the 
WA3 data. The original WA3 paper |13] provides the 
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results of the partial- wave analysis of the Kinr system in 
the form of plots for the intensity in the {K* {892)Tr)s, 
pK, KS{U30)tt, /o(1370)X, and {K*{892)Tr)D chan- 
nels, together with the phases of the corresponding am- 
plitudes, measured relative to the {K*{892)Tr)s ampli- 
tude. The uK data were not analyzed. In order to 
obtain an estimate of the parameters that enter the K- 
matrix model, we perform a fit of this model to the 
< < 0.05 GeV^ WAS data for the intensity of the 
m — Ktttt channels and the relative phases. Here \t'\ 
is the four momentum transfer squared with respect to 
the recoiling proton in the reaction K~p — > K~tt'^'k~p, 
and m denotes the magnetic substate of the Kmr sys- 
tem. Since the results of the analysis performed by the 
ACCMOR Collaboration are not sensitive to the choice 
of the value for the constant r in Eq. (O, we set r = 0. 
We seek solutions corresponding to positive values of the 
7± parameters, as found in the ACCMOR analysis [l^ . 
The data sample consists of 215 bins. The results of this 
fit are displayed in Fig. [1] and show a good qualitative 
agreement with the results obtained by the ACCMOR 
Collaboration [l^. We obtain = 855, with 26 free 
parameters, while the ACCMOR Collaboration obtained 
= 529. Although neither fit is formally a good one, 
the model succeeds in reproducing the relevant features 
of the data. 



C. Model for Ki production in B decays 

We apply the above formalism to the parameterization 
of the signal component for the production of Ki reso- 
nances in B decays. The propagation of the uncertain- 
ties in the if —matrix description of the ACCMOR data 
to the model for Ki production in B decays is a source 
of systematic uncertainty and is taken into account as 
described in Sec. I VII 

In order to parameterize the signal component for 
the analysis of B decays, we set the background am- 
plitudes D, whose contribution should be small in the 
non-diffractive case, to 0. The backgrounds arising from 
resonant and non-resonant B decays to the {K'K'k)'k fi- 
nal state are taken into account by separate components 
in the fit, as described in Sec. |Vl The parameters of K 
and the offset phases 6i are assumed to be independent 
from the production process and are fixed to the values 
extracted from the fit to WAS data (Table|T|. Finally, we 
express the production couplings fpa and fpb in terms of 
two real production parameters ^ = (t?, </>): fpa = cos-d, 
fpb = siniJe*'^, where e [0,7r/2], </> e [0,27r]. In this 
parameterization, tani^ represents the magnitude of the 
production constant for the i<ri(1270) resonance relative 
to that for the Ki{lAQQ) resonance, while is the relative 
phase. 

The dependence of the selection efficiencies and of the 
distribution of the discriminating variables (described in 
Sec.|V| on the production parameters C are derived from 
Monte Carlo studies. For given values of C,^ signal MC 




M (GeV) M (GeV) 

FIG. 1: Results of the fit to the < < 0.05 GeV^ WAS 
data. Intensity (left) and phase relative to the (T^"* (892)7r)s 
amplitude (right) for the (a) (Tf* (892)7r)s, (b, c) pK, (d, e) 
A'o(1430)7r, (f, g) /o(1370)A', and (h, i) (A'*(892)7r)D chan- 
nels. The points represent the data, the solid lines the total 
fit function, and the dashed lines the contribution from the 
background. 



samples for B decays to the [K-k-k)-!: final states are 
generated by weighting the {Kt:t:)7t population accord- 
ing to the amplitude ^^^^^^^^^(^^^N)^*' where the term 
{Kmrli) consists of a factor describing the angular dis- 
tribution of the Ktttt system resulting from Ki decay, 
an amplitude for the resonant tttt and Ktt systems, and 
isospin factors, and is calculated using the formalism de- 
scribed in Refs. p^, i2j]. The ujK channel is excluded 
from the sum, since the uj tt^tt^ branching fraction 
is only (1.5St°;i:^) 

compared to the branching fraction 
(89.2 ± 0.7) % of the dominant decay uj tt+tt'tt^ [i3]. 
Most of the Ki ojK decays therefore result in a dif- 
ferent final state than the simulated one. We account 
for the Ki — > ujK transitions with a correction to the 
overall efficiency. In Fig. [5] we show the reference frame 
chosen to evaluate the distributions of the products of 
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Parameter 


Fitted value 


Ma 


1.40 ± 0.02 


Mi 


1.16 ±0.02 


e 


72° ± 3° 


7+ 


0.75 ±0.03 


7- 


0.44 ± 0.03 


/a3 


0.02 ± 0.03 


fb3 


0.32 ±0.01 


fai 


-0.08 ±0.02 


fhi 


0.16 ±0.01 


fab 


0.06 ±0.01 


f 

Jb5 


U.zi zh U.U4 


52 


-31° ± 1° 


53 


82° ± 2° 


54 


78° ± 4° 


^5 


20° ± 9° 



spectively, in the Ki rest frame (Fig. [^b)- We define j3 
as the polar angle of the flight direction of m relative to 
the direction of the momentum of I (Fig. [^t). The re- 
sulting angular parts of the transition amplitudes for S-, 
P-, and D-wave decays of the Ki axial vector ( = l"*") 
mesons with scalar (J^ = 0+) and vector (J^ = 1^) 
intermediate resonances Xs_ci are given by: 



/ 3 

As — \ — (cos 6 cos /? + sin 9 sin (3 cos $) (15) 

V OTT 



(16) 



A-D = \ (-2 cos 9 cos /3 ± sin9sin/3cos$) . (17) 
V IBtt 



TABLE I: Parameters for the Tf-matrix model used in the 
analysis of B decays. 




Z P 



FIG. 2: Definition of (a) the coordinate axes in the Ki rest 
frame, (b) the angles O and $ in the Ki rest frame, and (c) the 
angle l3 in the rest frame of the Xs,d intermediate resonance. 



For the tttt and Kn resonances, the following ^-wave 
Breit-Wigner parameterization is used [2^ : 



BWim) = (tt) 



[mor(TO)]^/^ 



with 



(mg - 


- m?) 


— imQT{m) 


q{m) 


21+1 




_q{mo)_ 




1 + R^q^{m) 



(18) 



r(m) = r(mo)— 

m 

, (19) 
where toq is the nominal mass of the resonance, T{m) is 
the mass-dependent width, r(mo) is the nominal width 
of the resonance, q is the breakup momentum of the res- 
onance into the two-particle final state, and = 25 
GeV"2. The ir^(1430) and /o(1370) amplitudes are 
also parameterized as Breit-Wigner functions. For the 
/'irQ(1430) we assume a mass of 1.250 GeV and a width of 
0.600 GeV 119], while for the /o(1370) we use a mass of 
1.256 GeV and a width of 0.400 GeV ;.25]. This param- 
eterization is varied in Sec. IVII and a systematic uncer- 
tainty evaluated. 



B KiTT decays, where Ki decays proceed through the 
intermediate resonances Xs = {/C* (892), (1430)} or 
Xd = {Pi /o(1370), w}. Final state particles arc labeled 
with a subscript {k,l,m,n}, according to the following 
scheme: B° Kti^k, ^ X^^t:+ ^ K+tt- or 
B° ^ K+TT-, K+ X°K+, X° TT+TT- for neutral B 

meson decays, and 5+ -^i^fe j -^i ^ -^s^'"'/ ' -^t ^ 
Kln+ or B+ ^ i^? ^ X^^K^ X^ 7r+^^- for 

charged B meson decays. The angular distribution for 
the Ki system produced in B decays can be expressed 
in terms of three independent angles (9, /3, <&). In the 
Ki rest frame, we define the Y axis as the normal to 
the decay plane of the Ki^ and orient the Z axis along 
the momentum of / (Fig. [5^). 9 and $ are then the 
polar and azimuthal angles of the momentum of fc, re- 



IV. EVENT RECONSTRUCTION AND 
SELECTION 

The B^ —f K^n^ candidates are reconstructed in the 
K^TT^TT" decay mode by means of a vertex fit 
of all combinations of four charged tracks having a zero 
net charge. Similarly we reconstruct B^ K'^tt'^ can- 
didates, with K^TT+TT- ", by combining Kg candi- 
dates with three charged tracks. We require the recon- 
structed mass mK-n-n to lie in the range [1.1,1.8] GeV. 
Charged particles arc identified as either pions or kaons, 
and must not be consistent with the electron, muon or 
proton hypotheses. The candidates are reconstructed 
from pairs of oppositely-charged pions with an invariant 
mass in the range [486, 510] MeV, whose decay vertex is 
required to be displaced from the Ki vertex by at least 
3 standard deviations. 
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The reconstructed B candidates are characterized 
by two almost uncorrelated variables, the energy- 
substituted mass 



mEs = ■^(-s + po-Pb)VSo'-p| (20) 
and the energy difference 

AE = E*s- IV~S, (21) 

where (£'o,Po) and {Eb,Pb) are the laboratory four- 
momenta of the T{AS) and the B candidate, respectively, 
and the asterisk denotes the CM frame. We require 
5.25 < TTiEs < 5.29 GeV and \AE\ < 0.15 GeV. For 
correctly reconstructed B candidates, the distribution of 
mEs peaks at the B-meson mass and AE at zero. 

Background events arise primarily from random com- 
binations of particles in continuum e^e^ qq events 
(g = u, d, s, c). We also consider cross feed from other B 
meson decay modes than those in the signal. 

To separate continuum from BB events we use vari- 
ables that characterize the event shape. We define the 
angle 9t between the thrust axis [1^ of the B candidate 
in the T(4S') frame and that of the charged tracks and 
neutral calorimeter clusters in the rest of the event. The 
distribution of |cos6't| is sharply peaked near 1 for qq 
jet pairs and nearly uniform for _B-meson decays. We re- 
quire I cos6't| < 0.8. We construct a Fisher discriminant 
J- from a linear combination of four topological variables: 
the monomials Lq = T^iPi and L2 = X^i Pi I cos 6*,* ^ , 
I cos 6*^,1 and |cos6'|;| [23 . Here, p* and 0* are the CM 
momentum and the angle of the remaining tracks and 
clusters in the event with respect to the B candidate 
thrust axis. 9q and 9^ are the CM polar angles of the B- 
candidate thrust axis and i3-momentum vector, respec- 
tively, relative to the beam axis. In order to improve the 
accuracy in the determination of the event shape vari- 
ables, we require a minimum of 5 tracks in each event. 

Background from B decays to final states containing 
charm or charmonium mesons is suppressed by means 
of vetos. A signal candidate is rejected if it shares 
at least one track with a B candidate reconstructed in 
the B° D-n+, B° D*-tt+, B+ ^ D°t:+, or 
B+ D*°Tr~^ decay modes, where the D meson in the 
final states decays hadronically. A signal candidate is 
also discarded if any tt'^tt^ combination consisting of the 
primary pion from the B decay together with an oppo- 
sitely charged pion from the Ki decay has an invariant 
mass consistent with the cc mesons Xcoi^P) or Xci(l-P) 
decaying to a pair of oppositely charged pious, or J/ip 
and ip{2S) decaying to muons where the muons are mis- 
identified as pious. 

We define H as the cosine of the angle between the di- 
rection of the primary pion from the B decay and the nor- 
mal to the plane defined by the Ki daughter momenta in 
the Ki rest frame. We require |7i| < 0.95 to reduce back- 
ground from B — > Vn decay modes, where is a vector 
meson decaying to Kmr, such as _R'*(1410) or if* (1680). 



The average number of candidates in events containing 
at least one candidate is 1.2. In events with multiple 
candidates, we select the candidate with the highest 
probability of the B vertex fit. 

We classify the events according to the invariant 
masses of the tt+tt^ and K^tt" {K^tt^) systems in the 
{Ki) decay for B^ {B~^) candidates: events that sat- 
isfy the requirement 0.846 < rriKTr < 0.946 GeV belong 
to class 1 ( "if* band" ) ; events not included in class 1 for 
which 0.500 < m^^ < 0.800 GeV belong to class 2 ("p 
band"); all other events are rejected. The fractions of se- 
lected signal events in class 1 and class 2 range from 33% 
to 73% and from 16% to 49%, respectively, depending on 
the production parameters C- About 11% to 19% of the 
signal events are rejected at this stage. For combinato- 
rial background, the fractions of selected events in class 
1 and class 2 are 22% and 39%, respectively, while 39% 
of the events are rejected. 

The signal reconstruction and selection efficiencies de- 
pend on the production parameters C- For B° modes 
these efficiencies range from 5 to 12% and from 3 to 8% 
for events in class 1 and class 2, respectively. For i3+ 
modes the corresponding values are 4-9% and 2-7%. 

V. MAXIMUM LIKELIHOOD FIT 

We use an unbinned, extended maximum-likelihood 
(ML) fit to extract the event yields ns^r and the parame- 
ters of the probability density function (PDF) 'Ps,r- The 
subscript r = {1,2} corresponds to one of the resonance 
band classes defined in Sec. IIVI The index s represents 
the event categories used in our fit model. For the anal- 
ysis of i3° modes, these are 

1. signal, 

2. combinatorial background, 

3. 5° if*(1410)+7r-, 

4. 5° if*(892)"7r+7r- + p^K+n-, 

5. ^ ai(1260)±7r=F, and 

6. ^ D7.n+. 
For i?+ modes, these are 

1. signal, 

2. combinatorial background, 

3. B+ ^ if*(1410)"7r+, 

4. B+ -> if*(892)+7r+7r- +p^K'^tt+, and 

5. B+ -> K*{892)+p". 

The likelihood £e,r for a candidate e to belong to class 
r is defined as 

£e,r — ^ ^s,r ^s^rip^e: C; (^2) 
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where the PDFs are formed using the set of observables 
Xe = {Ai?, TTiES, ^ n^KTTTr, I'^l} and the dependence 
on the production parameters C is relevant only for the 
signal PDF. ^ represents all other PDF parameters. 

In the definition of £e,r the yields of the signal cate- 
gory for the two classes are expressed as a function of the 
signal branching fraction B as ni.i = B x Ng^ x ei(^) 
and ni^2 = B x Ngg x e2{C)j where the total selection 
efficiency er(C) includes the daughter branching fractions 
and the reconstruction efficiency obtained from MC sam- 
ples as a function of the production parameters. 

For the i?° modes we perform a negative log-likelihood 
scan with respect to d and 4>. Although the events in class 
r = 2 are characterized by a smaller signal-to-background 
ratio with respect to the events in class r — 1, MC studies 
show that including these events in the fit for the B'^ 
modes helps to resolve ambiguities in the determination 
of (j) in cases where a signal is observed. At each point of 
the scan, a simultaneous fit to the event classes r = 1, 2 
is performed. 

For the 5+ modes, simulations show that, due to a less 
favorable signal-to-background ratio and increased back- 
ground from B decays, we are not sensitive to 4> over a 
wide range of possible values of the signal BF. We there- 
fore assume (j> = tt and restrict the scan to -d. At each 
point of the scan, we perform a fit to the events in class 
r = 1 only. The choice (jj = tt minimizes the variations 
in the fit results associated with differences between the 
mKTTTT PDFs for different values of (j). This source of 
systematic uncertainty is accounted for as described in 
Sec. I VII The variations in the efficiency ei as a function 
of (j) for a given i) can be as large as 30 %, and are taken 
into account in deriving the branching fraction results as 
discussed in Sec. I VIII 

The fitted samples consist of 23167 events (S" modes, 
class 1), 38005 events {B^ modes, class 2), and 9630 
events modes, class 1). 

The signal branching fractions are free parameters in 
the fit. The yields for event categories s = 5, 6 (-6° 
modes) and s — 5 (5+ modes) are fixed to the values 
estimated from MC simulated data and based on their 
previously measured branching fractions [H, HI]. The 
yields for the other background components are deter- 
mined from the fit. The PDF parameters for combinato- 
rial background are left free to vary in the fit, while those 
for the other event categories are fixed to the values ex- 
tracted from MC samples. 

The signal and background PDFs are constructed as 
products of PDFs describing the distribution of each ob- 
servable. The assumption of negligible correlations in 
the selected data samples among the discriminating vari- 
ables has been tested with MC samples. The PDFs for 
AE and toes of the categories 1, 3, 4, and 5 are each pa- 
rameterized as a sum of a Gaussian function to describe 
the core of each distribution, plus an empirical function 
determined from MC simulated data to account for the 
tails of each distribution. For the combinatorial back- 
ground we use a first degree Chebyshev polynomial for 



AE and an empirical phase-space function (29| for toes: 



f{x) cx x^/l - x2 exp - x"^)] , (23) 

where x = 2toes/'\/s and is a parameter that is de- 
termined from the fit. The combinatorial background 
PDF is found to describe well both the dominant quark- 
antiquark background and the background from random 
combinations of B tracks. 

For all categories the J- distribution is well described 
by a Gaussian function with different widths to the left 
and right of the mean. A second Gaussian function with 
a larger width accounts for a small tail in the distribution 
and prevents the background probability from becoming 
too small in the signal region. 

The rriK-K-K distribution for signal depends on C,. To 
each point of the C scan, we therefore associate a differ- 
ent nonparametric template, modeled upon signal MC 
samples reweighted according to the corresponding val- 
ues of the production parameters and (p. Production 
of K*{141Q) and ai(1260) resonances occurs in B back- 
ground and is taken into account in the rriKjnr and \Ti.\ 
PDFs. For all components, the PDFs for |7i| are param- 
eterized with polynomials. 

We use large control samples to verify the tties, AE, 
and PDF shapes, which are initially determined from 
MC samples. We use the i?" — > D~Tr'^ decay with 
D- ^ K+TT-TT- , and the B+ D°tt- decay with 1?° 
A'gTT+TT", which have similar topology to the signal 5° 
and 5+ modes, respectively. We select these samples by 
applying loose requirements on uies and AE, and requir- 
ing for the D candidate mass 1848 < to^,- < 1890 MeV 
and 1843 < mjja < 1885 MeV. The selection require- 
ments on the B and D daughters are very similar to 
those of our signal modes. These selection criteria are 
applied both to the data and to the MC events. There 
is good agreement between data and MC samples: the 
deviations in the means of the distributions are about 
0.5 MeV for niES, 3 MeV for AE, and neghgible for JF. 



VI. SYSTEMATIC UNCERTAINTIES 

The main sources of systematic uncertainties are sum- 
marized in Table |TT1 For the branching fractions, the 
errors that affect the result only through efficiencies are 
called "multiplicative" and given in percentage. All other 
errors are labeled "additive" and expressed in units of 
10^6. 

We repeat the fit by varying the PDF parameters ^, 
within their uncertainties, that are not left floating in the 
fit. The signal PDF model excludes fake combinations 
originating from misreconstructed signal events. Poten- 
tial biases due to the presence of fake combinations, or 
other imperfections in the signal PDF model, are esti- 
mated with MC simulated data. We also account for 
possible bias introduced by the finite resolution of the 
{■d, (j)) likelihood scan. A systematic error is evaluated 
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by varying the _K'i(1270) and ii'i(1400) mass poles and 
ii'— matrix parameters in the signal model, the param- 
eterization of the intermediate resonances in Ki decay, 
and the offset phases 5i. We test the stability of the fit 
results against variations in the selection of the "X*" 
and "p bands," and evaluate a corresponding system- 
atic error. An additional systematic uncertainty origi- 
nates from potential peaking BB background, including 
B -> i-Ta* (1430)77 and B ii:*(1680)7r, and is evaluated 
by introducing the corresponding components in the def- 
inition of the likelihood and repeating the fit with their 
yields fixed to values estimated from the available ex- 
perimental information [l^. We vary the yields of the 
B° ^ ai(1260)±7r=F and 5° ^ D'+^-^-T^^ (for the B" 
modes) and 5+ K*~^p^ (for the B'^ modes) event 
categories by their uncertainties and take the resulting 
change in results as a systematic error. For i?+ modes, 
we introduce an additional systematic uncertainty to ac- 
count for the variations of the parameter. The above 
systematic uncertainties do not scale with the event yield 
and are included in the calculation of the significance of 
the result. 

We estimate the systematic uncertainty due to the in- 
terference between the B Kitt and the B K*TnT + 
pKn decays using simulated samples in which the decay 
amplitudes are generated according to the results of the 
likelihood scans. The overall phases and relative con- 
tribution for the K*t:t: and pKir interfering states are 
assumed to be constant across phase space and varied 
between zero and a maximum value using uniform prior 
distributions. We calculate the systematic uncertainty 
from the RMS variation of the average signal branching 
fraction and parameters. This uncertainty is assumed 
to scale as the square root of the signal branching frac- 
tion and does not affect the significance. The system- 
atic uncertainties in efficiencies include those associated 
with track finding, particle identification and, for the _B+ 
modes, Kg reconstruction. Other systematic effects arise 
from event selection criteria, such as track multiplicity 
and thrust angle, and the number of B mesons. 



TABLE II: Estimates of systematic errors, evaluated at the 
absolute minimum of each — ln£ scan. For the branching 
fraction, the errors labeled (A) , for additive, are given in units 
of 10~^, while those labeled (M), for multiplicative, are given 
in percentage. 





B 


K 


+ - 




Quantity 


B 


■& 


4> 


B 


■d 


PDF parameters (A) 


0.8 


0.01 


0.15 


1.4 


0.07 


MC/data correction (A) 


0.8 


0.00 


0.01 


1.0 


0.02 


ML fit bias (A) 


0.6 


0.03 


0.02 


2.0 


0.08 


Fixed phase (A) 


— 


— 


— 


0.6 


0.06 


Scan (A) 


0.9 


0.04 


0.16 


0.0 


0.04 


Ki ii— matrix parameters (A) 


2.2 


0.01 


0.36 


0.5 


0.05 


Ki onset phases (A) 


0.2 


0.01 


0.02 


0.0 


0.00 


Ki intermediate resonances (A) 


0.5 


0.00 


0.06 


0.2 


0.02 


K* /p bands (A) 


0.2 


0.05 


0.00 


1.2 


0.05 


Peaking BB bkg (A) 


0.8 


0.01 


0.13 


1.0 


0.01 


Fixed background yields (A) 


0.0 


0.00 


0.00 


0.4 


0.02 


Interference (A) 


6.0 


0.25 


0.52 


10.6 


0.43 


MC statistics (M) 


1.0 






1.0 




Particle identification (M) 


2.9 






3.1 




Track finding (M) 


1.0 






0.8 




Kg reconstruction (M) 








1.6 




cos 6't(M) 


1.0 






1.0 




Track multiplicity (M) 


1.0 






1.0 




Number BB pairs (M) 


1.1 






1.1 





Fig. [5ji as a function of the production parameters for 
B^ and i?+ modes, respectively. The associated statisti- 
cal error <jb{C,) ^f each point given by the change in 
B when the quantity —2 ln£(S; C,) increases by one unit, 
is displayed in Fig. [5^ and Fig. [Sf. Systematics are in- 
cluded by convolving the experimental two-dimensional 
likelihood for d and (/),£ = £(Sniax(C); C); with a two- 
dimensional Gaussian that accounts for the systematic 
uncertainties. In Fig. [B^ and Fig. we show the result- 
ing distributions in and (j). The 68% and 90% prob- 



VII. FIT RESULTS 

Figures [3] and U] show the distributions of AE, toes 
and mKiT-n for the signal and combinatorial background 
events, respectively, obtained by the event-weighting 
technique s'Plot [3^. For each event, signal and back- 
ground weights are derived according to the results of 
the fit to all variables and the probability distributions 
in the restricted set of variables in which the projection 
variable is omitted. Using these weights, the data are 
then plotted as a function of the projection variable. 

The results of the likelihood scans are shown in Ta- 
ble mil and Fig. [5j At each point of the C scan the 
— 21n£(S;C) function is minimized with respect to the 
signal branching fraction B. Contours for the value 
^max(C) tha-t maximizes C{B\ C,) are shown in Fig.[SJ: and 



TABLE III: Results of the ML fit at the absolute minimum 
of the — In £ scan. The first two rows report the values of the 
production parameters (j)) that maximize the likelihood. 
The third and fourth rows are the reconstruction efficiencies, 
including the daughter branching fractions, for class 1 and 
class 2 events. The fifth row is the correction for the fit bias 
to the signal branching fraction. The sixth row reports the re- 
sults for the B 7^1(1270)^ + A'i(1400)7r branching fraction 
and its error (statistical only). 







B+ ^ KU+ 




0.86 


0.71 




1.26 


3.14 (fixed) 


ei (%) 


3.74 


1.36 


£2 (%) 


1.68 




Fit bias correction (xlO~®) 


4-0.0 


+0.7 


B (xlO-®) 


32.1 ±2.4 


22.8 ± 5.1 
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rngs (GeV) A E (GeV) m^^^ (GeV) 



FIG. 3: sPlot projections of signal onto rriES (left), A_E (center), and mK-mr (right) for class 1 (top), class 2 (middle), 
and B^ class 1 (bottom) events: the points show the sums of the signal weights obtained from on-resonance data. For tties 
and AiJ the solid line is the signal fit function. For mK-mr the solid line is the sum of the fit functions of the decay modes 
iCi(1270)7r + _ft'i(1400)7r (dashed), if*(1410)7r (dash-dotted), and Jf* (892)7r7r (dotted), and the points are obtained without 
using information about resonances in the fit, i.e., we use only the rriES, Ai5, and T variables. 




niEs (GeV) A E (GeV) m^^^ (GeV) 

FIG. 4: sPlot projections of combinatorial background onto rriES (left), A_E (center), and mjf^^ (right) for B'^ class 1 (top), 
class 2 (middle), and B^ class 1 (bottom) events: the points show the sums of the combinatorial background weights obtained 
from on-resonance data. The solid line is the combinatorial background fit function. For mK-mr the points are obtained without 
using information about resonances in the fit, i.e., we use only the rriEs, AE, and T variables. 



14 




FIG. 5: (a, b) — ln£ scan (systematics not included) in the production parameters i9 and for the (a) B'^ and (b) modes. 
The cross in (a) indicates the position of the absolute minimum in the — ln£ scan. A second, local minimum is indicated 
by a star and corresponds to an increase in A(— ln£) of 2.7 with respect to the absolute minimum, (c, d) Contours for the 
B /ri(1270)7r + i(:i(1400)7r branching fraction (in units of 10"^) extracted from the ML fit for the (c) B° and (d) B+ modes, 
(e, f) Contours for the statistical error (in units of 10~^) on the B — » /S'i(1270)7r + ifi(1400)7r branching fraction for the (e) 
B° and (f) B+ modes. 
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ability regions are shown in dark and light shading, re- 
spectively, and are defined as the regions consisting of all 
the points that satisfy the condition £(r) > x, where the 
value X is such that ]c.{r)>x '^(^' 4>)'^'&^4> = 68% (90%). 
The significance is calculated from a likelihood ratio test 
A(— 21n£) evaluated at the value of i? that maximizes 
the likelihood averaged over (j>. Here A(— 21n£) is the 
difference between the value of — 21n£ (convolved with 
systematic uncertainties) for zero signal and the value 
at its minimum for given values of ^. We calculate the 
significance from a distribution for A(— 21n£) with 
2 degrees of freedom. We observe nonzero K^tt~ 
and K^tt~^ branching fractions with 7.5(t and 3.2cr 

significance, respectively. 

We derive probability distributions for the B 
ifi(1270)7r + Ki{UOO)t:, B i(:i(1270)7r, B 

ii'i(1400)7r, and B Kiatt branching fractions. 

At each point in the C plane we calculate the distri- 
butions for the branching fractions, given by f{B;C) = 
cC{B;^), where c is a normalization constant. Sys- 
tematics are included by convolving the experimental 
one-dimensional likelihood C{B; C) with a Gaussian that 
represents systematic uncertainties. Branching frac- 
tion results are obtained by means of a weighted av- 
erage of the branching fraction distributions defined 
above, with weights calculated from the experimental 
two-dimensional likelihood for and (fi. 

For each point of the ^ scan the B — > i4ri(1270)7r, 
B — > iiri(1400)7r, and B — > Kiatt branching fractions 
are obtained by applying ^-dependent correction factors 
to the B ifi(1270)7r + ifi(1400)7r branching frac- 
tion associated with that ^ point. The correction fac- 
tor is calculated by reweighting the signal MC samples 
by setting the production parameters {fpa, fpb) equal to 
(0, e**^ sini?), (cosi?, 0), and (l/p^l cos0, — l/p^l sin6'), for 
B ifi(1270)7r, B ifi(1400)7r, and B Kiatt, re- 
spectively, where fpA = cos cos 9 — e*"^ sin -d sin 6 and 9 
is the Ki mixing angle [l^, for which we use the value 
9 = 72° (see Table U. 

From the resulting distributions f{B) we calculate 
the corresponding two-sided intervals at 68% probabil- 
ity, which consist of all the points B > that sat- 
isfy the condition f{B) > x, where x is such that 
IfiB)>x B>o/(^)^^ = 6^^"- The upper limits (UL) at 
90% probability are calculated as /o<g<(7i f{B)dB — 
90%. The results are summarized in Table IIVI (statis- 
tical only) and Table IVl (including systematics). 
We measure 6(5° ^ Xi(1270)+7r--|-Xi(1400)- 



q 1 -1-0. o 
■^■^-0.7 



and BiB+ 

-1+2.9 



Xi(1270)07r-i 



10- 

Ki{U00fT:+) = 2.9li;^ x 10"^ (< 8.2 x 10"^), where 
the two-sided ranges and upper limits are evaluated at 
68% and 90% probability, respectively, and include sys- 
tematic uncertainties. 

Including systematic uncertainties we obtain 
two-sided intervals (in units of 10^"''): B{B'^ 
ifi(1270)+7r-) e [0.6,2.5], B{B° Ki(1400)+7r- 



the 



[0.8,2.4], B{B° 



e [0.4,2.3], B{B^ 



i4:i(1270)°7r+) e [0.0,2.1] (< 4.0), BiB+ 
is:i(1400)°7r+) e [0.0,2.5) (< 3.9), B{B+ K°^tt+) G 
[0.0,2.1] (< 3.6), where the two-sided ranges and the 
upper limits are evaluated at 68% and 90% probability, 
respectively. 



VIII. BOUNDS ON |Aa| 

We use the measurements presented in this work to 
derive bounds on the model uncertainty \Aa\ on the 
weak phase a extracted in B'-* ai(1260)*7rT de- 
cays. We use the previously measured branching frac- 
tions of B° ai(1260)±7r=F, B° ai(1260)-i4:+ and 
B+ ai{1260)+K° decays 4] and the CP-violation 
asymmetries ,6] as input to the method of Ref . fs"] . The 
values used are summarized in Tables I VII and IVIII 

The bounds are calculated as the average of |Aq;|+ = 



off 



-a I and jAaj 



c^off ~ which are obtained from 



the inversion of the relations 



cos2(a^ — a) > 



l~A 



cos2(Q;^ff - a) > 



2i?+ 



± 2 
CP 

+ 
± 



1 - A 



± 2 
CP 



(24) 



where we have defined the following ratios of CP- 
averaged rates 



Rt 



R 



+ — 



-\yiBiK%7:+: 
fl^B{aU-) 

}?flB{atK^) 
flB{a^^+) ■ 



The CP asymmetries A^p in B^ 



afiT^ decays are 
related to the time- and flavor-integrated charge asym- 
metry A^'^ i by 



A 



CP 



AC) + C 



AllpC 



A(jp — 



AC) 



AC 

~C 



Ac'pC - 



AC 



C and AC parameterize the flavor-dependent direct CP 
violation and the asymmetry between the CP-averaged 
rates B{aiTr~) and B{a^n'^), respectively Q: 



C± AC 



1^^ 
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TABLE IV: Branching fraction results for B — > Knr decays, in units of 10~®, and corresponding confidence levels (C.L., 
statistical uncertainties only). For each branching fraction we provide the mean of the probability distribution, the most 
probable value (MPV), the two-sided interval at 68% probability, and the upper limit at 90% probability. 



Channel 


Mean 


MPV 


68% C.L. interval 


90% C.L. UL 


B" 7fi(1270) + 7r- + 7fi(1400) + 7r~ 


3.2 


3.1 


(2.9,3.4) 


3.5 


^ A'i(1270) + 7r- 


1.7 


1.6 


(1.3,2.0) 


2.1 


B° A:i(1400) + 7r- 


1.6 


1.6 


(1.3,1.9) 


2.0 




1.5 


1.4 


(1.0,1.9) 


2.2 


B+ /fi(1270)"7r+ + is:i(1400)%+ 
B+ 7^1(1270)%+ 
B+ ^-1(1400)%+ 
B+ Kl^^+ 


2.9 
1.1 
1.8 
1.1 


2.3 
0.3 
1.7 
0.2 


(1.6,3.5) 
(0.0, 1.4) 
(1.0,2.5) 
(0.0,1.5) 


4.5 

2.5 
2.0 
2.3 



TABLE V: Branching fraction results for B — » /sTitt decays, in units of 10~^, and corresponding confidence levels (C.L., 
systematic uncertainties included). For each branching fraction we provide the mean of the probability distribution, the most 
probable value (MPV), the two-sided interval at 68% probability, and the upper limit at 90% probability. 



Channel 


Mean 


MPV 


68% C.L. interval 


90% C.L. UL 


B" K-i{l27i)) + 'K- + Ki(1400) + 7r- 


3.3 


3.1 


(2.4,3.9) 


4.3 


^ A'i(1270) + 7r~ 


1.7 


1.7 


(0.6,2.5) 


3.0 


B° ^ Ji:i(1400) + 7r- 


1.6 


1.7 


(0.8,2.4) 


2.7 


B° ^ i^+^TT- 


1.6 


1.4 


(0.4, 2.3) 


2.9 


B+ ->■ JTi (1270)%+ + if 1 (1400)%+ 
B+ 7^1 (1270)%+ 
B+ ^-1(1400)%+ 
B+ if^^+ 


4.6 
1.7 
2.0 
1.6 


2.9 
0.0 
1.6 
0.2 


(1.2,5.8) 
(0.0,2.1) 
(0.0,2.5) 
(0.0,2.1) 


8.2 
4.0 
3.9 
3.6 
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TABLE VI: Summary of the branching fractions used as input 
to the calculation of the bounds on |Aq| i4]. 



Decay 


mode 




Branching fraction 
(in units of 10"^) 




ai(1260)± 




33.2 ±3.8 ±3.0 




ai(1260)- 


K+ 


16.3 ± 2.9 ± 2.3 


B+ ^ 


ai(1260)-* 




33.2 ±5.0 ±4.4 



TABLE VII: Summary of the values of the CP— violation pa- 
rameters used as input to the calculation of the bounds on 
|Aa| i. 



Quantity 


Value 




-0.07 ±0.07 ±0.02 


S 


0.37 ±0.21 ±0.07 


AS 


-0.14 ±0.21 ±0.06 


C 


-0.10 ±0.15 ±0.09 


AC 


0.26 ±0.15 ±0.07 



where the decay amplitudes for B^{B ) ^ ai(1260)^7r^ 
are 

A+ = A{B° a+TT") , EE A{B" a^7r+) , 

A+ = A(^ ar7r+) , A- = A(B° a+Tr") . 

The CP-averaged rates are calculated as 

Bia+n-) = i6(a±^T)(l ± AC ± A^c'pC), 

B{a-7:+) = i6(a±7rT)(l - AC - All^C), 

where B{a^'K^) is the flavor-averaged branching fraction 
of neutral B decays to ai(1260)=^7rT 

For the constant A = |Ks|/|Kd| = | V'cd | / 1 | we 
take the value 0.23 [13] ■ The decay constants Jk — 
155.5 ± 0.9 MeV and = 130.4 ± 0.2 MeV O are 
experimentally known with small uncertainties. For the 
decay constants of the oi and Kia mesons the values 
fa, = 203 ± 18 MeV ^ and Jk.^ = 207 MeV [3 are 
used. For /kia we assume an uncertainty of 20 MeV. 
The value assumed for the Jkia decay constant is based 
on a mixing angle 9 = 58° Q, because /kia is not avail- 
able for the value 9 = 72° used here (see Table H]); this 
discrepancy is likely accommodated within the accuracy 
of the present experimental constraints on the mixing 
angle. Using naive arguments based on SU(3) relations 
and the mixing formulae, we have verified that the de- 
pendence of fxiA on the mixing angle is rather mild in 
the 6 range [58,72]°. It should be noted that due to a 
different choice of notation, a positive mixing angle in 
the formalism used by the ACCMOR Collaboration [l^ 
and in this paper corresponds to a negative mixing angle 
with the notation of Ref . @ . 



We use a Monte Carlo technique to estimate a prob- 
abihty region for the bound on \aes — a\. All the CP- 
averaged rates and CP-violation parameters participat- 
ing in the estimation of the bound are generated accord- 
ing to the experimental distributions, taking into account 
the statistical correlations among A^p, C, and AC [1^. 

For each set of generated values we solve the system 
of inequalities in Eq. (j24p . which involve jo;^ — a\ and 
\a~g — a\, and calculate the bound on jofcff ~ ct\ from 

]acff -a\< (Ja+f - a] ± \a-g ~ al)/2. (25) 

The probability regions are obtained by a counting 
method: we estimate the fraction of experiments with 
a value of the bound on joioff — ot\ greater than a given 
value. We obtain \a^s - a\ < 11°(13°) at 68% (90%) 
probability. 

The determination of aoff Q presents an eightfold am- 
biguity in the range [0°,180°]. The eight solutions are 
aeff = (11 ± 7)°, aeff = (41 ± 7)°, aeff - (49 ± 7)°, 
aeff = (79 ± 7)°, aeff = (101 ± 7)°, = (131 ± 7)°, 
a^s = (139 ± 7)°, aoff = (169 ±7)° @. Assuming that 
the relative strong phase between the relevant tree am- 
plitudes is negligible [1] it is possible to reduce this am- 
biguity to a twofold ambiguity in the range [0°,180°]: 
aoff = (11 ± 7)°, aeff = (79 ± 7)°. We combine the so- 
lution near 90°, aoff = (79 ± 7)° [f|, with the bounds on 
jaoff — a] and estimate the weak phase a = (79±7±11)°. 
This solution is consistent with the current average value 
of a, based on the analysis of P ^ tttt, B pp, and 
B ^ pn decays [Hi^. 

IX. SUMMARY 

We present results from a branching fraction measure- 
ment of P ~> i4:i(1270)7r and A'i(1400)7r decays, ob- 
tained from a data sample of 454 million T(4S') PP 
events. The signal is modeled with a if -matrix formal- 
ism, which accounts for the effects of interference be- 
tween the isri(1270) and i4ri(1400) mesons. Including 
systematic and model uncertainties, we measure B{B^ — s- 
i4:i(1270)+7r- ± Xi(1400)+7r-) = i-lt^j x 10^^ and 
B{B+ Xi(1270)°7r+±ifi(1400)°7r+) = 2.9ti;? x 10"^ 
(< 8.2 X 10"^ at 90% probabihty). A combined signal for 
the decays P° ifi(1270)+7r- and P° ifi(1400)+7r- 
is observed with a significance of 7.5ct, and the fol- 
lowing branching fractions are derived for neutral P 
meson decays: 6(P° -> Xi(1270)+7r-) e [0.6,2.5] x 
10-^ i3(P" Xi(1400)+7r-) e [0.8,2.4] x 10"^, and 
6(P° K+^n-) e [0.4,2.3] x 10~^, where the two- 
sided intervals are evaluated at 68% probability. A sig- 
nificance of 3.2cr is obtained for P+ fs:i(1270)°7r+ ± 
ifi(1400)°7r+, and we derive the following two-sided in- 
tervals at 68% probability and upper limits at 90% proba- 
bility: B{B+ /i:i(1270)°7r+) € [0.0, 2.1] x 10"^ (< 4.0x 
10-^), B{B+ ifi(1400)07r+) e [0.0,2.5) x 10"^ (< 
3.9 X 10-5), and B{B+ K^^7r+) e [0.0,2.1] x 10"^ (< 
3.6 X 10-5). 
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Finally, we combine the results presented in this paper 
with existing experimental information to derive an in- 
dependent estimate for the CKM angle a, based on the 
time-dependent analysis of CP-violating asymmetries in 
B° ai(1260)±7r=F, and find a = (79± 7± 11)°. 

X. ACKNOWLEDGEMENTS 

We thank Ian Aitchison for helpful discussions and 
suggestions. We are grateful for the extraordinary con- 
tributions of our PEP-II colleagues in achieving the 
excellent luminosity and machine conditions that have 
made this work possible. The success of this project 
also relies critically on the expertise and dedication of 
the computing organizations that support BaBAR. The 
collaborating institutions wish to thank SLAC for its 



support and the kind hospitality extended to them. 
This work is supported by the US Department of En- 
ergy and National Science Foundation, the Natural Sci- 
ences and Engineering Research Council (Canada), the 
Commissariat a I'Energie Atomique and Institut Na- 
tional de Physique Nucleaire et de Physique des Partic- 
ules (France), the Bundesministerium fiir Bildung und 
Forschung and Deutsche Forschungsgemeinschaft (Ger- 
many), the Istituto Nazionale di Fisica Nucleare (Italy), 
the Foundation for Fundamental Research on Matter 
(The Netherlands) , the Research Council of Norway, the 
Ministry of Education and Science of the Russian Fed- 
eration, Ministerio de Educacion y Ciencia (Spain), and 
the Science and Technology Facilities Council (United 
Kingdom). Individuals have received support from the 
Marie-Curie lEF program (European Union) and the A. 
P. Sloan Foundation. 



[1] V. Laporta, G. NarduUi and T. N. Pham, Phys. Rev. 

D 74, 054035 (2006); Erratum ibid. D 76, 079903 (2007). 
[2] G. Calderon, J. H. Munoz and C. E. Vera, Phys. Rev. 

D 76, 094019 (2007). 
[3] H.-Y. Cheng and K.-C. Yang, Phys. Rev. D 76, 114020 

(2007) . 

[4] BABAR CoUaboration, B. Aubert et al, Phys. Rev. Lett. 

97, 051802 (2006); Phys. Rev. Lett. 100, 051803 (2008). 
[5] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 

99, 261801 (2007); Phys. Rev. Lett. 99, 241803 (2007). 
[6] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 

98, 181803 (2007). 

[7] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); 
M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 
652 (1973). 

[8] M. Gronau and J. Zupan, Phys. Rev. D 73, 057502 

(2006); Phys. Rev. D 70, 074031 (2004). 
[9] M. Gronau and J. L. Rosner, Phys. Lett. B 595, 339 

(2004). 

[10] M. Beneke, M. Gronau, J. Rohrer, and M. Spranger, 

Phys. Lett. B 638, 68 (2006). 
[11] BABm Collaboration, B. Aubert et al, Phys. Rev. D 76, 

052007 (2007). 

[12] Particle Data Group, C. Amsler et al, Phys. Lett. B 667, 
1 (2008). 

[13] Except as noted explicitly, we use a particle name to 
denote either member of a charge conjugate pair. 

[14] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. 
B 254, 288 (1991). 

[15] Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 
87, 161601 (2001); Belle Collaboration, H. Yang et al, 
Phys. Rev. Lett. 94, 111802 (2005); BABm Collabora- 
tion, B. Aubert et al, Phys. Rev. Lett. 101, 161801 

(2008) . 

[16] BABMi Collaboration, B. Aubert et al, Nucl. lustrum. 

Methods A479, 1 (2002). 
[17] G. Benelli et al. Nuclear Science Symposium Conference 



Record, 2005 IEEE, 2, 1145 (2005). 

[18] S. Agostinelli et al, Nucl. lustrum. Methods Phys. Res., 
Sect. A 506, 250 (2003). 

[19] ACCMOR Collaboration, C. Daum et al, Nucl. 
Phys. B187, 1 (1981). 

[20] I. J. R. Aitchison, Nucl. Phys. A189, 417 (1972). 

[21] M. Nauenberg and A. Pais, Phys. Rev.126, 360 (1961). 

[22] ACCMOR Collaboration, C. Daum et al, Nucl. 
Phys. B182, 269 (1981). 

[23] M. G. Bowler, M. A. V. Game, I. J. R. Aitchison, and J.B. 
Dainton, Nucl. Phys. B97, 227 (1975); ACCMOR Collab- 
oration, C. Daum et al, Nucl. Phys. B182, 269 (1981); 
M. G. Bowler, J. Phys. G3, 775 (1977); L. Stodolsky, 
Phys. Rev. Lett. 18, 973 (1967); M. G. Bowler, J. Phys. 
G3, 1503 (1977). 

[24] D. J. Herndon, P. Soding, and R. J. Cashmore, Phys. 
Rev. D 11, 3165 (1975). 

[25] Particle Data Group, R. L. Kelly et al. Rev. Mod. Phys. 
52, SI (1980). 

[26] A. de Riijula, J. Ellis, E. G. Floratos, and M. K. Gaillard, 

Nucl. Phys. B 138, 387 (1978). 
[27] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 70, 

032006 (2006). 

[28] Heavy Flavor Averaging Group, E. Barberio 

et al, [arXiv:0808.1297' and online update at 

http: / /www.slac. stanford.edu/xorg/hfag 
[29] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. 

B 241, 278 (1990). 
[30] M. Pivk and F. R. Le Diberder, Nucl. lustrum. Methods 

Phys. Res., Sect. A 555, 356 (2005). 
[31] J. C. R. Bloch, Yu. L. Kalinovsky, C. D. Roberts, and 

S. M. Schmidt, Phys. Rev. D 60, 111502 (1999). 
[32] J. Charles et al (CKMFitter Group), Eur. Phys. Jour. 

C 41, 1 (2005); M. Bona et al (UTfit Collaboration), 

JHEP 0603, 080 (2006). 



